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We find the Gutzwiller projected Fermi sea wave function(GWF) has the correct phase structure 
to describe the Icink nature of the doped hoies in the ground state of the one dimensional t — J 
model. We find the failure of the GWF for general value of J/t and electron density n can be 
attributed to the residual charge correlation in the ground state. We find such residual charge 
correlation is well described by a XXZ-type effective Hamiltonian. Based on these observations, a 
Pfaffian-type variational wave function is proposed and is found to reproduce correctly the global 
phase diagram and corresponding correlation functions of the one dimensional t— J model, including 
the Luther-Emery phase in the low electron density and large J/t region. 

PACS numbers: 74.20.Mn,74.25.Ha,75.20.Hr 



I. INTRODUCTION 

The Gutzwiller projected wave functions(GWFs) are 
widely used to approximate the ground state of the t— J 
model and the Heisenberg model. In these models, lo- 
cal electronic correlation, as manifested in the no double 
occupancy constraint of electrons, plays a vital role in 
determining the low energy physics. Such strong local 
correlation make these systems difficult to study analyti- 
cally. In the variational approach based on GWFs, these 
models are first treated in the mean field approximation 
in which the local constraint is relaxed to a global one. 
Then the local constraint is enforced afterwards by the 
Gutzwiller projection which simply filters out the un- 
physical components with doubly occupied sites in the 
mean field state. 

The above variational strategy is used extensively in 
the study of the high temperature superconductors and 
quantum antiferromagnets. After many years of efforts, 
it is now believed that the Gutzwiller projected d-wave 
BCS state describe well the superconducting state of the 
high temperature superconductors ^1, ^ ,^ ^]. Quite re- 
cently, progress is also made on the understanding of the 
quasiparticle properties above such a state[5|,|a,l7J,|3| . The 
same kind of wave function is also used in recent studies 
on the exotic orders and exotic excitations of frustrated 
quantum antiferromagnet|9l ll0l UlL Il2 | . 

An unresolved issue about the GWF is that it is not 
clear if the posteriorly executed projection can capture 
the kinematic effect of the local constraints, even quali- 
tatively. In this paper, we address this issue with the one 
dimensional t — J model. 

The one dimensional t~ J model has been studied ex- 
tensively by a broad band of methods including Bethe- 
Ansatz solution, conformal field thcoriesjlj, \^A lla|j 

quantum Monte carlo llfil . exact diagonalization llTl . and 
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l22l l23l |24| . Many properties concerning the ground state 
of this model are now well established. This give us the 
unique opportunity to judge the validity of a given ap- 
proximation. The one dimensional i — J model is ex- 
actly soluble at J/t = and J/t = 2. For J/t = 0, 
the spin and the charge degree of freedom of the system 
are totally separatedjlg- The spin part is described by 
the Heisenberg model on the squeezed chain with doped 
holes removed, while the charge part is described by a 
noninteracting spinless Fermion system. For J/t = 2, the 
system is supersymmetric and it is found that the GWF 
provides a fairly accurate approximation for the ground 
state of the svstem jl,4l22j . For general value of J/t and 
electron density n, the system is a Tomonaga-Luttinger 
liquid(TLL) below a critical value Jc/t around 2.5. The 
correlation exponent of the TLL varies continuesly with 
J/t and njrgl- For J/t > Jc/t, the system is unstable 
toward phase separation. For small n and J/t > 2, there 
is also a small region in which the system exhibit a spin 
gapmH. 

The Gutzwiller projected Fermi sea wave function and 
its variants has long been used to describe the ground 
state of the one dimensional t~ J model. It is well known 
that this wave function provides an excellent description 
of the undoped case of the model, namely the spin ^ 
Heisenberg spin chain[ig|. However, the same wave func- 
tion is not that satisfactory for the doped system, except 
for the supersymmetric case of J/t = 2. For example, it 
fails to predict the TLL behavior in the small J/t region. 
A 2]<ip peak in the spin structure factor is also missed by 
this wave function. Since the wave function is parameter 
free, it also gives no clue on the origin of the spin gap 
state and the phase separation at large J/t. 

It is generally believed that the problems with the 
GWF originate from the insufficient account of the charge 
correlation in the system. Along this line of think- 
ing, various kind of Jastrow factor are proposed to rem- 



edy the drawbacks of GWF. For example, Hellberg and 
Mele introduced a long range Jastrow factor of the form 
\F{ri^,rji)\'' and succeed in reproducing the TLL be- 
havior, where \F(rir,ri\)\'^ is a Slater determinant of all 
the electron positions |20|. Yokoyama and Ogata found a 
short range repulsive Jastrow factor is able to restore the 
2ki? peak in the spin structure factor, while a sufficiently 
attractive Jastrow factor can cause phase separation [2^. 
However, both wave functions have difficulties in repro- 
ducing the correct phase diagram. For example, the spin 
gap state is missed in both wave functions. At the same 
time, both wave functions predict a fully phase-separated 
state along the boundary of phase separation, which is 
in fact an oversimplification|l7j . More importantly, no 
understanding on the physical origin of the proposed Jas- 
trow factor is available and it is hard to judge if a similar 
modification is relevant for higher dimensional system. 

For the sake of possible extension to higher dimensional 
case, it is important to know the reason that the simple 
GWF fails before any modification on it is made. As 
mentioned above, it is the residual charge correlation in 
the system which is responsible for the failure of GWF. 
In this paper we make this statement more precise by 
showing that the GWF has the correct phase structure 
to describe the kink nature of the doped holes in the 
ground state of the one dimensional t ~ J model. In 
fact, we find the spin structure factor of the GWF in the 
squeezed chain coordinate is almost identical to that of a 
half filled spin chain. Thus the missing 2ki? peak in the 
spin structure factor for small J /t should be recovered if 
the removed holes are reinserted into the squeezed chain 
in the right manner. 

The physical origin of the residual charge correlation 
can be easily seen if one reformulate the GWF in terms of 
the slave Boson theory [23, 123 • In the slave Boson theory, 
the constrained electron operator is decomposed &s c[ ^ = 

fi^bi, in which fj^ is a spin ^ neutral Fermion called 
spinon and bi is a spinless charge 1 Boson called holon. 
The local constraint now takes the form of an equality, 
X^cr /i CT/i.cr+^l^j ~ ^- In tcrms of thc slavc Boson thcory, 
the GWF corresponds to a state with all holon condensed 
into the zero momentum state. However, the holon is 
not a true Boson but a hard core Boson as a result of 
the local constraint. For general value of J/t, there is 
also an effective attraction between the holons caused by 
the exchange term of the t — J model. Thus a XXZ-type 
effective Hamiltonian should be a good approximation 
for the residual charge correlation. 

Based on these observations, a Pfaffian-type varia- 
tional wave function is proposed for the ground state 
of the one dimensional t ~ J model. This wave func- 
tion, which has only one parameter, reproduce well the 
global phase diagram of the model, including the Luther- 
Emery(LE) phase in the small n and large J/t region. It 
is found that this wave function also reproduces well var- 



ious correlation functions of the system and provides a 
refined picture for the phase separation at large J/t. 

The paper is organized as follows. Section II is devoted 
to the investigation of the properties of the GWF. In 
Section III, the new variational scheme and the Pfaffian- 
type wave function are introduced. The phase diagram 
and correlation functions determined from this new vari- 
ational wave function are presented in Section IV. Sec- 
tion V summarize the paper and includes a discussion on 
related issues. 



II. 



THE GWF 



The one dimensional t — J model reads 
n^-t ^(cj^c^+i,^ + h.c.) -f J ^(Si • Si+i - -1 



HiH+l)i 



(1) 

in which S.^ = i Y^ap cla'^apCip and Ui = J2a ^la^ia- The 
electron in this model is subjected to the constraint of no 
double occupancy 



Xl^Lcia < 1. 



(2) 



The ground state of the one dimensional t — J model is 
governed by a well defined phase structure. This can be 
most easily seen at half filling when the system reduces to 
the Heisenberg spin chain. For the Heisenberg model, it 
is well known that the ground state satisfy the Marshall 
sign rule|27l l2a |. The rule says that the ground state 
wave function is real in the Ising basis and its sign is 
given by (— 1)^^ up to a global phase, where iV| denotes 
the number of down spins in the even sublattice. This 
sign rule is a manifestation of the antiferromagnetic spin 
correlation in the ground state. With such a sign rule, 
one easily verify that (S^ • Sj) < for i and j belonging to 
different sublattices and (Si-Sj) > for i and j belonging 
to the same sublattice. 

The ground state at finite doping is governed by a sim- 
ilar sign rule. It can be easily checked that all matrix 
elements of the t — J Hamiltonian are negative definite 
in a wave function that satisfy the Marshall sign rule on 
the squeezed chain. The squeezed chain is the chain in 
which the sites occupied by the doped holes are removed. 
This can be seen by noting that the motion of holes in 
this model do not disturb the spin configuration on the 
squeezed chain. Thus, the ground state of the one dimen- 
sional t—J model should satisfy such a modified Marshall 
sign rule. With such a modified Marshall sign rule, one 
easily see that the holes in the ground state behaves as 
an antiphase domain wall for spin. 

Now we show that the GWF satisfies the Marshall sign 
rule on the squeezed chain. The GWF reads 



|GWF)=n(l-"'T"a)|FS), 



(3) 



in which |FS) denotes the simple Fermi sea. In the nat- 
ural basis Y\- ■ cLcL |0), the amplitude of GWF is given 
by the following Vandermont determinant 



Q</3 



l<7n 



in which ^Tw = cxp[— iki?(^^ «„ + J2iJi)] is ^ plane 
wave factor, Zi^ = exp(i^^^) and Zj, = exp{i-j^) are 
chord coordinates of up spins and down spins. 

Now we exchange the up spin at site ii and the down 
spin at site ji . The resultant change in phase is given by 



A$ 



arg 



(n 



a>l 



z, 



n^ 



z,. 



1>1 



Zi, — z. 



(5) 



ji 



Since \Z\ = 1, arg(^ 



-) is nothing but the angle in 



the segment Zi-^ — Zj^ in the unit circle. Noting the fact 
that in a circle the angles in the same segment equal one 
another and sum of the opposite angles of quadrilaterals 
equals TT, one easily find that A$ = A^'cTt, in which N^ de- 
notes the number of electrons between site i\ and site j\. 
Taking into account the sign due to Fermion exchange, 
one find the change in phase is in accordance with the 
modified Marshall sign rule. Following essentially the 
same steps, one can also verify the case of exchanging a 
hole and an electron. 



approximation for J jt = 2, it is natural to expect that 
the same behavior to hold for arbitrary J/t and n. 

Two conclusions follow directly from the above rea- 
soning. First, since the spin correlation on the squeezed 
chain is already well described by the GWF, the missing 
2kp peak in the spin structure factor should be recov- 
ered if the removed holes are correctly reinserted into 
the squeezed chain, or, the the missing 2kp peak should 
be attributed to the residual charge correlation in the 
system. Second, since the squeezed spin chain picture 
is argued to hold for arbitrary J/t and n, a single wave 
function may suffice to describe the whole phase diagram 
of the one dimensional t — J model, including the spin 
gap phase at small n and large J/t. 



III. 



THE NEW VARIATIONAL SCHEME 



The origin of the residual charge correlation can be 
most easily seen by reformulating the GWF in terms of 
the slave Boson theory. In the slave Boson theory, the 
constrained electron operator is decomposed as c] ^ — 

fi^bi, in which fj^ represents the Fermionic spinon and 
bi represents the Bosonic holon. In terms of these slave 
particles, the t — J model reads 

Ti. = Ti-t + "Hj 



squeezed chain 
Heisenberg 




FIG. 1: Spin structure factor of the GWF in the squeezed 
coordinate as compared with that of a half filled spin chain. 
The calculation is done on a 204 sites lattice which is quarter 
filled. 

Thus the GWF has the right phase structure to de- 
scribe the ground state of the one dimensional i— J model 
and the kink nature of the doped holes in it. In fact, this 
conclusion can be made even stronger. In Figure 1, we 
plot the spin structure factor of the GWF in the squeezed 
coordinate and compare it with that of a half filled spin 
chain. We see the two are almost identical with each 
other. Since the spin degree of freedom is described ex- 
actly by a Heisenberg model on the squeezed chain at 
J/t = 0, while the GWF provides an exceedingly good 



^* = -tT.^fln+i.^blib7 + h.c.) 

■i,cr 
i 

in which S/ = i Ea/3 //a^«/3/»/3 and nj = Ea //«/»«■ 
The no double occupancy constraint now takes the form 
of an equality 



E/l/-+^^^' = i- 



(6) 



When the local constraints Eq.(6) is exactly satisfied, the 
factor biblbi+ibl_^i appearing in Tij plays no role and can 
be neglected. 

In the mean field treatment, an RVB order parameter 
X = Eq < fi+iafia > IS introduccd to decompose the 
interaction term. At the same time, the local constraint 
is relaxed to a global one. The mean field Hamiltonian 
for the spinon and the holon part read 26J 

nf = -(tx + ^)J2iflf,+,. + h.c.) 
n' = -txJ2{blb,+i + h.c), 



in which x is the hole density. The mean field ground 
state is given by the product of the spinon Fermi sea and 



the holon Bose condensate 



1$) 



ibl 



'k=oJ 



n /kt/kMo). 



(7) 



k<kj7 



When this state is projected into the subspace that sat- 
isfy the constraint Eq.(6), we get the GWF. 

In the mean field theory, the holon is a free Boson and 
condenses in the ground state. However, due to the local 
constraint, the holon is actually a hard core Boson which 
can not condense in one spatial dimension. The uncon- 
densed nature of the hard core Boson in Id originates 
from the kinematic effect of the local constraint. Due 
to this constraint, the Hilbert space for the one dimen- 
sional hard core Boson system becomes disconnected at 
the single particle level. We note for comparison that the 
Hilbert space of the spinon part is still connected even 
when the local constraint is enforced. Thus the holon 
should be treated as hard core Boson rather than free 
Boson. 

Another source of the residual charge correlation is 
provided by the superexchange term of the t — J model. 
When two electrons are next to each other, they enjoy 
an attraction due to the superexchange. This attraction 
is not captured by the mean field order parameter x a-nd 
should be reintroduced. 

Combining these considerations, the residual charge 
correlation beyond the GWF should be described by the 
following XXZ-type effective Hamiltonian 



W„ = - J2^blk+i + h.c.) -vJ2 im+A+ik 



(8) 



^t 



in which 6] is the operator for hard core Boson and v is 
the rescaled attraction. If we denote the ground state of 
Hy as A„, then PgAi,|GWF) should be a good variational 
wave function for the one dimensional t — J model. 

Although Tiv is exactly soluble [22|, an explicit form 
for A^ is available only in limited cases. For v — 0, 
Ay is nothing but the Hellberg-Mele Jastrow factor with 
1^=1. For w = 1, A^ is a constant and our proposed wave 
function reduce to the GWF. At the quarter filling, Hv 
exhibit particle-hole symmetry. In a separated paper |2fl|, 
we show a Hellberg-Mele-type variational wave function 
provides an exceedingly good description for the ground 
state of the XXZ model in the S"^ = sector. How- 
ever, away from the particle-hole symmetric point, the 
Hellberg-Mele wave function cease to be a good approx- 
imation. 

For general value of v and Boson density, we have to 
resort to approximation. Through the Jordan- Wigner 
transformation, the XXZ Hamiltonian can be rewritten 



Ti-v 



E 



(c|cj+i -I- h.c.) 



.\^JJ 






^i^i+l^i+'^'^ii 



(9) 



in which c\ is a spinless Fermion. For this Hamiltonian, 
we adopt the BCS approximation to decouple the interac- 
tion term. The BCS ground state for the spinless Fermion 
reads 



nK+«k4c-k)io)' 



(10) 



k>0 



in which ^^ = -^, Ak = Asin(k), Ck = -2cos(k) - /^ 
and Ek = yeJTA^. Here A is the BCS gap for the 
spinless Fermion and is treated as the only variational 
parameter in our theory(the chemical potential /i can be 
determined by the density equation and is not an inde- 
pendent parameter). In real space, the BCS state for the 
spinless Fermion takes the form of a Pfaffian. A Pfaffian 
is a square root of the determinant of a antisymmetric 
matrix of even order[31j. In our case, the matrix element 
of the antisymmetric matrix is given by 



fzj =^— sin(k(i-j)), 
k>o"i^ 



(11) 



in which i and j denote the coordinates of the spinless 
Fermions. Thus our variational wave function for the one 
dimensional t — J model is given by 



* = Pf(A)|GWF), 



(12) 



in which Pf(A) is the Pfafhan for the holes which are 
now spinless Fermions. 



IV. 



RESULTS 



A. Ground state phase diagram 

The ground state phase diagram determined from the 
Pfaffian- type wave function is presented in Figure 2. 

The phase diagram contains three distinct phases. For 
small and intermediate value of J/t, the system is in the 
TLL phase in which both charge and spin excitation are 
gapless. For larger value of J/t, the system is unsta- 
ble towards phase separation. At small n and large J/t, 
there is a small region in which the system exhibits a spin 
gap. In the spin gap phase, the charge excitation is still 
gapless. Following the convention, this phase is termed 
Luther-Emery liquid. 

The phase boundaries are determined as follows. To 
illustrate the idea, we plot the variational energy as a 
function of the electron density for J/t — 2 and J/t — 3 
in Figure 3. For J/t — 2, the energy curve is concave 
everywhere so that a homogenous phase is globally stable 
for all electron density. For J/t = 3, a convex region 
appears at intermediate values of electron density in the 
energy curve. In this case, the boundaries for the globally 
stable phases are given by the two tangency points shown 
in the figure, while the boundaries for the locally stable 
phases are given by the two inflexion points. 
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FIG. 2: Ground state phase diagram of the one dimensional 
t ~ J model determined from the Pfafhan-type variational 
wave function. The dotted lines indicate the boundaries for 
the existence of locally stable phases. Here TLL denotes 
Tomonaga-Luttinger liquid, LEL denotes Luther-Emery liq- 
uid, while PS denotes phase separated state. 
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points, the system is unstable towards phase separation. 
The density of the phase separated phases are given by 
two tangency points. For 2.5 < J/t < 3.2, the system 
phase separates into a hole rich phase and an electron 
rich phase. For 3.2 < J/t < 3.4, the hole rich phase 
is replaced by a empty phase. For J/t > 3.4, a fully 
phase separated state is realized in which the electron 
rich phase is replaced by a half filled spin chain. 

The convex region of the energy curve diminishes to 
zero at about J/t = 2.5. The phase boundary between 
the TLL phase and the LEL phase for J/t < 2.5 is de- 
termined by examining the infrared behavior of the spin 
structure factor S{q). In the spin gap phase, S{q) should 
be quadratic at small g, while in the TLL phase a lin- 
ear behavior is exDected|32l|. For the charge excitation, 
a similar criteria exists on the density structure factor 
Niq). 

The existence of the LEL phase is quite unexpected 
from the point of view of the mean field theory. In the 
mean field theory, the spinon is still described by a filled 
Fermi sea which is by definition gapless. However, after 
Gutzwiller projection the spinon get entangled with the 
holon. Such an entanglement change drastically the spin 
correlation of the system. 



B. Correlation functions 

Four correlation functions are evaluated in this work. 
They are the momentum distribution function defined as 



(^)""?^5I^^*>^J'^)^ 



ikin-Tj) 



(13) 




the spin structure factor defined as 



FIG. 3: Variational energy per site e as a function of the 
electron density Ue for J/t = 2 and J/t — 3. For clarity's sake, 
a linear decreasing background of the energy is subtracted. 
a — ^ji^lne^o is the initial slope of the energy curve. The 
arrows above the curve indicate the locations of the inflexion 
points, while the arrows below the curve indicate the locations 
of the tangency points. The phase boundaries are determined 
from these points as explained in the text. 



For electron density that lies between the two tangency 



Sik) = jTT.(^!S])e^''^^'' 



r,) 



«J 



the charge structure factor defined as 



«j 



and the pair correlation function defined as 



(14) 



(15) 



(16) 



in which A^ is the annihilation operator for a nearest- 
neighboring pair 

Aj = ^(ci|Ci+ix - ciia+i]). (17) 

First we present the result for the TLL phase. The 
correlation functions for J/t = 0, 1 and 2 at quarter fill- 
ing are shown in Figure 4. For comparison's sake, we 



also plot the result calculated from the Hellberg-Mele 
wave function. From the figure we see that the corre- 
lation functions calculated from the Pfaffian-type wave 
function are almost identical with that calculated from 
the Hellberg-Mele wave function, apart from the small 
deviations due to critical fluctuations. Since the Pfaffian 
is derived from a BCS mean field approximation in which 
a gap opens up, the residual charge correlation described 
by it is short ranged. Thus the Pfaffian-type wave func- 
tion should exhibit Fermi-liquid behavior, as is clear in 
Figure 4. To recover the critical fluctuations, one should 
go beyond the mean field approximation. 






FIG. 5: (a) The momentum distribution function n(fc), (b) 
the spin structure factor S{k), (c) the charge structure factor 
C{k), and (d) the singlet pairing correlation function P{k) 
at quarter filling for J/t=2.5(black square) and J/t=3(red up 
triangle). The solid lines denote the result calculated from 
the Hellberg-Mele variational wave function. 



wave function is slightly lower than that of the Hellberg- 
Mele wave function. For larger value of J/t, the ordering 
is reversed. However, both wave functions give good esti- 
mate for the ground state energy before phase separation. 



FIG. 4: (a) The momentum distribution function n(fc), (b) 
the spin structure factor S{k), (c) the charge structure factor 
C(fc), and (d) the singlet pairing correlation function P{k) at 
quarter filling for J/t=0(black square), J/t=l(red up trian- 
gle), and J/t=2(green circle). The solid lines denote the result 
calculated from the Hellberg-Mele variational wave function. 

For J/t > 2, the Pfafhan-type wave function becomes 
less satisfactory for the quarter filled system. In Figure 5, 
we plot the correlation functions of the quarter filled sys- 
tem at J/t = 2.5 and 3, the latter of which is very close to 
the boundary of phase separation. Near the boundary of 
phase separation, the Hellberg-Mele wave function starts 
to develop charge instability, as is clear from Fig. 5(c). 
This tendency is missed by the Pfaffian-type wave func- 
tion. Instead, a structure at 2kp remains evident in the 
correlation functions. This is to be expected, since we 
start from a Fermionic description of the residual charge 
correlation. In fact, it is quite amazing that the BCS 
approximation remains to be a good approximation for 
J/t as high as 2.5(the optimized value for the BCS gap 
A is approximately given by J/t at quarter filling). 

To quantify the above discussion, we plot in Figure 6 
the relative error in the variational energy for both the 
Pfaffian-type wave function and the Hellberg-Mele wave 
function. For small J/t, the energy of the Pfaffian-type 
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FIG. 6: Relative error in variational energy at quarter filling. 
The exact value of the ground state energy is taken from ,2^ . 

Although the Hellberg-Mele wave function provides a 
good description for the quarter filled system, it fails 
badly at low electron density. On the other hand, 
the Pfafhan-type wave function describes quite well the 
physics in the low density regime, including the spin gap 



phase at large J/i. To illustrate this, we plot in Figure 7 
the error in variational binding energy for a single pair of 
electrons calculated from both wave functions. From the 
figure we see the Pfaffian-type wave function is almost 
exact for all values of J /t in the low density limit. We 
think this explains why the spin gap phase can be cor- 
rectly reproduced by the Pfaffian-type variational wave 
function. 




j/t 



FIG. 7: Error in variational binding energy for a single pair 
of electrons calculated from the Pfaffian-type wave function 
and the Hellberg-Mele wave function. The inset shows an 
expanded view of the < J/t < 2 region. 



Now we present the correlation functions for the LEL 
phase at small n and large J/t. In Figure 8, the corre- 
lation functions for J/t = 2.8 and n = 0.06, a system 
deeply inside the LEL phase, are plotted. As mentioned 
above, the spin gap manifests itself in the quadratic be- 
havior of the spin structure factor in the small q limit. We 
note that spin gap state is metastable in a much larger 
region than that of the LEL phase. 



V. SUMMARY AND DISCUSSION 

In this paper, we have carried out a variational study 
of the one dimensional t — J model. We find the failure 
of the simple GWF should be attributed to the resid- 
ual charge correlation. Reformulating the GWF in terms 
of the slave Boson theory, we find the residual charge 
correlation should be described by a XXZ-type effective 
Hamiltonian. Based on this observations, a Pfaffian-type 
variational wave function is proposed for the one dimen- 
sional t — J model. We find this wave function, which has 
only one variational parameter, reproduces correctly the 
global phase diagram and the corresponding correlation 
functions. 




FIG. 8; Correlation functions for J/t — 2.8 and n — 0.06, 
a system deeply inside the spin gap phase, (a) The momen- 
tum distribution function n{k), (b) the spin structure factor 
S{k), (c) the charge structure factor C{k), and (d) the singlet 
pairing correlation function P{k). 



It is interesting to note the way in which the spin cor- 
relation is affected by the charge degree of freedom in this 
model. Through the investigation of the phase structure 
of the ground state wave function, we find the doped holes 
behave as anti-phase domain walls for the spin correla- 
tion. We show further the spin degree of freedom of the 
system is well approximated by a half filled spin chain in 
the squeezed coordinates throughout the phase diagram. 
For small electron density, the effect of the charge degree 
of freedom on the spin part can be so drastic as to in- 
duce a gap in the excitation spectrum of the latter. This 
spin gap is beyond mean field description and should be 
attributed to the strong entanglement of spin and charge 
degree of freedom in the projected subspacc. 

It is also interesting to note the effect of the local con- 
straint for this system. In conventional GWF, the effect 
of the local constraint is taken into account posteriorly 
by filtering out the unphysical components in the unpro- 
jected state. In this paper, we find this procedure may 
fail when the kinematic effect of the local constraint is 
essential for establishing (or, more accurately, destroying) 
the mean field correlation in the unprojected state. The 
one dimensional t—J model provides a particular exam- 
ple of this type. In the one dimensional t — J model, 
the Hilbert space for the charge degree of freedom is 
disconnected at the single particle level due to the lo- 
cal constraint. When the local constraint is relaxed, the 
connectivity of the Hilbert space for the charge degree 
of freedom is changed in a qualitative manner. Such a 
change in the connectivity of the Hilbert space is respon- 
sible for the appearance of the Bose condensation of the 
charged particle in the mead field theory and is ultimately 
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responsible for the failure of the GWF to describe the 
Tomonaga-Luttinger behavior of the system. 

For a full understanding of the residual charge corre- 
lation in the one dimensional t — J model, one should 
also take into account the attraction due to the exchange 
term. In the mean field treatment, the exchange term is 
decoupled in the fj^fja channel which can not account 
for such a charge correlation effect. We find this attrac- 
tion counteracts with the effect of the local constraint 
and cancels it out around J/t = 2. This explains the 
excellentness of the GWF at the supersymmetric point. 

The critical behavior of the one dimensional t—J model 
is not correctly described by the Pfaffian-type wave func- 
tion. This is natural since the Pfafhan itself is derived 
from a BCS-type mean field treatment. Although single 
particle condensation is gone in this treatment, a conden- 
sate of pair of spinless Fermion still exists. To recovery 
the correct critical behavior, one must get rid of such a 
condensate. One way to achieve this is to introduce a 
second charge correlator of the Hellberg-Mele-type, as is 
done in [2J| and [22|. According to our analysis, the prob- 
lem of finding a good variational description for the one 
dimensional t—J model reduces to that for the much sim- 
pler one dimensional XXZ model. We think the correct 
critical behavior should be recovered by a more accurate 
guess for the ground state wave function of the latter 
model. 

Finally, we mention possible generalization of the idea 
used in this work to the study of the two dimensional 
t — J model. In two spatial dimension, the kinematic 
effect of the local constraint should be less dramatic since 
the connectivity of the Hilbert space is not affected by 
the local constraint. This can also be seen from the fact 
that the two dimensional XXZ model does undergo Bose 
condensation at zero temperature. However, it is much 
subtler to analyze the interplay between the spin and 
charge degree of freedom in two spatial dimension since 
the two frustrate each other. Thus, the validity of the 
simple GWF in two dimension remains to be seen. 

This work is supported by NSFC Grant No.90303009. 
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